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In this paper we describe and compare difTerent methods used for accurate determination of forces 
acting on matter-wave packets in optical lattices. The quantum interference nature responsible for 
the production of both Bloch oscillations and coherent delocalization is investigated in detail. We 
study conditions for optimal detection of Bloch oscillation for a thermal ensemble of cold atoms 
with a large velocity spread. We report on the experimental observation of resonant tunneling in an 
amplitude-modulated (AM) optical lattice up to the sixth harmonic with Fourier-limited linewidth. 
We then explore the fundamental and technical phenomena which limit both the sensitivity and the 
final accuracy of the atomic force sensor at 10~^ precision level [J], with an analysis of the coherence 
time of the system and addressing few simple setup changes to go beyond the current accuracy. 

PACS numbers: 37.10.Jk, 03.75.-b, 03.75.Lm, 37.25. -^k, 04.80.-y 
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I. INTRODUCTION 

Coherent manipulation of cold atoms trapped in opti- 
cal lattice potentials finds powerful applications in dif- 
ferent fields like high resolution metrology quan- 
tum information processes Q, strongly correlated quan- 
tum phases as much as in the determination of 
forces at micrometer resolution In particular, 
loading cold atoms into tilted optical lattices led to 
the possibility to detect Bloch oscillations in the pres- 
ence of constant forces [l^ [T3. as well as Wannier- 
Stark ladders and tunneling while the introduc- 
tion of dynamical driving of the lattice depth has been 
used for several purposes: performing spectroscopy stud- 
ies of intcrband excitations [l6l - [l8j . characterizing the 
Mott-insulator regime [l^, quantum simulations of spin 
chains [2^ [2l| and demonstrating Loschmidt echoes in 
thermal gases [22l . [23l and their utilization as a high fi- 
delity atom mirror [23| . 

In the present paper we aim to show how it is pos- 
sible to merge the two phenomena of Bloch oscillations 
in tilted optical lattices and dynamical driving of lat- 
tice phase and amplitude to perform measurements of 
forces with high sensitivity and accuracy. This combina- 
tion leads to the realization of novel technique that we 
call delocalization-enhanced Bloch oscillations (DEBO), 
which is characterized in detail both theoretically and ex- 
perimentally in this work. We provide for the first time 
an unitary treatment of the dynamical system originated 
by an accelerated periodic potential which is driven by 
both amplitude (AM) or phase (PM) modulation. We 
give comprehensive theoretical explanation of the effect 
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originating from an interplay between external acceler- 
ation and the driving fields. This theoretical analysis 
and its predictions can be used to understand and inter- 
pret the experimental demonstrations about controlling 
and tailoring quantum transport of matter waves through 
coherent delocalization by resonant phase [l^l or ampli- 
tude [23} modulation, observation of quantum transport 
over macroscopic distances, macroscopic Bloch oscilla- 
tions in real space by off-resonant modulation (25| . and 
coherent phase imprint over the matter waves showing 
Loschmidt echoes. Furthermore, starting from the anal- 
ysis of space and time symmetries, we provide a demon- 
stration of the resonance spectrum accuracy of the co- 
herent delocalization technique. 

Particular attention is then paid to the utilization of 
these techniques to develop a force sensor for precision 
measurements of the gravitational acceleration g based 
on an accurate determination of Bloch frequency. We 
provide experimental measurements of Bloch frequency 
by direct observation of the modulation frequency at 
which resonant tunneling occurs, and we compare it with 
the one derived from DEBO technique in momentum 
space after several thousands periods. Since these two 
techniques rely on coherences of matter waves, we experi- 
mentally investigate the decoherence phenomena limiting 
the force measurement sensitivity. Finally we review all 
the important systematic effects and the specific impact 
on the two different techniques, providing the error bud- 
get and thus, for the first time, a significant comparison 
of their performance. 

The paper is organized as follows. We outline in Sec.lTTI 
the important aspects characterizing the quantum mo- 
tion in modulated optical lattices in terms of Wannier- 
Stark states, discussing it from the two alternative points 
of view of position and momentum space. We calculate 
the transport resonance spectrum by means of the Flo- 
quet theory [26[ . We provide an analytical simple model 
explaining the interference peak enhancement of Bloch 
oscillation visibility resulting from prior lattice modu- 
lation. We describe in Sec. IIIII the experimental setup 



2 



and procedures adopted in order to tune the governing 
parameters and study the system's observables. We ex- 
pUcitly discuss in Sec.|lV]the general advantages of using 
modulation driving to improve phase sensitivity in the 
measurement of Bloch oscillations frequency. The anal- 
ysis and comparison of sensitivity obtained with Bloch 
oscillation and resonant tunneling measurements is pre- 
sented in Sec. |Vl where the error budget is analyzed. 
Sec. IVII finallv contains concluding remarks and perspec- 
tives on how to further improve sensitivity for precision 
force measurements. 



II. THEORETICAL FRAMEWORK: QUANTUM 
MOTION IN MODULATED OPTICAL LATTICES 

A. The system: tilted optical lattices 

We consider the quantum motion of non-interacting 
ultracold atoms of mass M, loaded into the one- 
dimensional (ID) optical lattice U{z) 



C/(z) = --cos -z 



(1) 



with lattice constant d and lattice depth Uq [27|. This 
is originated by the interference pattern of two counter- 
propagating laser beams with wavelength = 27r/fc/^; 
the interference contrast controls directly the depth [/q, 
while the wavelength determines the lattice period d = 
Ai,/2. As we are interested in a non-dissipative potential, 
the laser beam frequency has to be far detuned from any 
atomic transition to make the scattering rate at which 
atoms spontaneously emit photons negligible. 

Motional degrees of freedom can be driven by exter- 
nal forces F{z) = —dzV{z) and quantum motion can be 
tailored by modulating the optical lattice either in am- 
plitude or in phase. The effective Hamiltonian governing 
the system is 

n{z,p, t) = ^ + + «/ W] - Pzfit) + V{z) , (2) 

where f{t) = sin(a;A/t — 0) is the time-dependent har- 
monic driving function, a and /3 are respectively the am- 
plitude and phase modulation depths. In the reference 
frame co-moving with the lattice (Kramers-Henneberger 
transformation), /? is a function of the spatial driving 
amplitude zq and ujm as /3 = M zoLoj^j. 

In the absence of both external forces and modula- 
tions, the stationary solutions of ^ arc Bloch states 
\ipB{k)) labelled by band index j and quasi-momcntum fc, 
which satisfy Bloch's theorem expressing invariance after 
translation of n lattice sites [l3| . The corresponding en- 
ergy spectrum consists of Bloch bands Ej{k). These are 
conveniently represented within the first Brillouin zone 
BZ = [— fc/,, fci], so that each band of index j is contin- 
uous over this interval. Within a semiclassical approach. 



the derivative of the energy band determines as a func- 
tion of k the group velocity of the wave packets in the 
optical lattice: 



Vg{k) = 



1 dEjjk) 
h dk 



(3) 



Depending on the particular quasi-momentum k distri- 
bution, an initially localized particle can propagate from 
site to site by means of the so-called Bloch tunneling, 
resulting in a ballistic spread over the whole lattice |28| . 
By applying external forces and lattice modulations, it 
is possible to control and coherently drive the quantum 
motion of wave packets. 

In the simplest case of a constant force F{z) = F, 
atoms undergo oscillations with angular frequency 



Fd 
IT 



(4) 



which manifest themselves both in momentum and real 
space; this oscillatory phenomenon is well known as 
Bloch oscillations (BO). The external constant force F 
breaks the translational symmetry, so that the Bloch 
states |'0s(fc)) are no longer eigenstates of H(a = 0,/3 = 
0). Nonetheless, we can still use them as a basis set. If in- 
terband Landau-Zener (LZ) tunneling can be neglected, 
one finds the time-evolved wavefunction 



{k~Ft/h)\i,i{k)) (5) 



where Gj{k) are periodic functions with period 27r/c?, 
which are determined by the initial state at i = 0. Thus, 
in case LZ tunneling can be neglected, the system be- 
haves as if the quasi-momcntum k evolves in time ac- 
cording to the semiclassical equation of motion 



dkjt) 
dt 



F. 



(6) 



Here, it is implicitly understood that the evolution of 
quasi-momentum is considered folded within the first 
Brillouin zone. 

The linear increase of the quasi-momentum results in 
a periodic change of the atomic momentum distribu- 
tion with the angular frequency when the quasi- 
momentum reaches the boundary of the Bloch band it is 
Bragg reflected [1^ . 

An alternative description of BOs can be given in terms 
of Wannier-Stark (WS) states l^'nj), which are quasi- 
stationary states constructed from the energy band of 
index j and centered on the lattice site of index n [soj . 
In essence, the —Fz potential breaks the translational 
symmetry, and the state energies at different sites become 
mismatched to produce a ladder 



f 



Ej — nFd ~ Ej — nhujB , 



(7) 



where the single energy step is determined by the Bloch 
frequency in Eq. (|4]), and Ej = d/ {2tt) J^^^dk Ej{k) rep- 
resents the average energy value of the j-th band. 
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Bloch tunneling is thus suppressed and the phe- 
nomenon of Wannier-Stark locahzation occurs, though 
non-zero probabihty for LZ tunnehng introduces a finite 
hfetime |3]| . 

If we consider a state with initial defined wave vector 
fco at t = in the lowest band, its time evolution will be 



n— — oo 

CO 



(8) 



where k{t) = ko + ujBt/d should be intended modulus 
the Brillouin zone. The latter result is identical to that of 
Eq. ([5]) previously derived from the semiclassical equation 
of motion. This second point of view puts in evidence the 
interferometric character of this phenomenon, where a 
phase imprint cogt occurs from site to site. Thus, observ- 
ing Bloch oscillations requires the initial state to be pre- 
pared with coherences between separated sites, namely 
in a superposition of at least two different WS states, 
e.g., I^'jij) and |^'„+i^j). Conversely, if the initial state is 
given by a single WS state, then the dynamics would sim- 
ply reduce to the accumulation of a trivial global phase. 
Site-to-site coherence implies in other terms that the co- 
herence length of the atomic cloud must be longer than 
the lattice period d, which for a thermal gas is of the 
order of the thermal de Broglie wavelength. 

Therefore, intuitively a way to enhance the visibility of 
BO momentum peak is to set a coherent broadening of 
the initial wave packet. We will come back to this point 
in Sec. Ill El where we show how amplitude and phase 
modulation of the optical lattices serve to the observation 
of DEBO's. 



B. Driving by nearly resonant amplitude and 
phase modulation 

We study the dynamics resulting from the time- 
dependent Hamiltonian in We show here that both 
the amplitude modulation (AM) (a ^ 0, /3 = 0) and 
phase modulation (PM) (a = 0, /3 7^ 0) of an optical 
lattice are very fiexible tools to steer quantum transport 
of atomic wave packets in tilted optical lattices, i.e., in 
the presence of an external constant force F. In par- 
ticular, we are interested in coherent delocalization of 
matter waves by means of resonant tunneling. Reso- 
nant tunneling occurs when the atoms absorb or emit 
energy quanta resonant with the ^-th harmonics 

of Bloch frequency ihujB, while they tunnel upwards or 
downwards along the optical lattice, respectively. For 
this reason we will consider nearly resonant modulation 
frequencies wa/, i.e. ujm — ^ws. Energy quanta are 
exchanged between the atoms and the photons of the 
modulated optical lattice, suggesting an interpretation 



of this transport mechanism in terms of photon-assisted 
tunneling [ll,!!!,!!!. 

In the rest of this work we will consider only physical 
regimes where interband LZ tunneling is fully negligible. 
For clarity, we therefore drop the band index j, since un- 
der this assumption the atomic dynamics remains con- 
fined in the same initial lattice band at all times. 

In the case of AM driving, we can use the WS basis to 
rewrite the Hamiltonian in Eq. ^ in the tight-binding 
form 



-f 00 



"Ham = -"^Wi3|*n>(*n| 



(9) 



+ CX3 



^ ( ^ C^M - 0) I Vl/„+,) (*„ I + h.C. 



with the coefficients Cf^ = (*„+£| cos(2fcLz)|^'„) 
representing the real-valued overlap integrals between 
resonantly-coupled WS states |5'n). In writing the 
Hamiltonian (|9]) we neglected all time-dependent cou- 
pling terms, except for those which are nearly reso- 
nant. For this reason we disregarded the interband 
coupling terms, too, as we assume the modulation fre- 
quency LdM be much smaller than the energy gap be- 
tween lattice bands. It is possible to express the Hamil- 
tonian in a time-independent form by transforming it 
into the rotating frame through the operator U^f ~ 
exp[-iY,n=-oo'^i^M/^)t\'ifn){'in\]- Usiug the rotating- 
wave approximation where the fast oscillating terms are 
neglected, one obtains 



-)-oo 



/ /tAM 



(10) 



where S = ujm — (-ujb is the detuning, the prime refer to 
the transformed reference frame, and the complex- valued 
coefficients represent the tunneling rates 



AM 



*„+^|cos(2fcLz)|*„) 



(11) 



Both parameters a and (j) can be precisely controlled in 
an experiment, offering great flexibility in steering atomic 
transport. For instance, a simple shift by tt of the phase 
(j) allows the tunneling rates J^^^ to be inverted |23| . 

It is worth noticing that the coefficient Cf^ can be 
viewed as the Franck-Condon factor of a stimulated two- 
photon Raman transition [34j . leading to the interpreta- 
tion of a two-photon Raman process underling the reso- 
nant tunneling among intraband WS states 135| . 

In the case of PM driving, following a procedure similar 
to the one used for AM modulation, one can express the 
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Hamiltonian in the tight-binding form 



State evolution 



n— — OO 

+ 00 

- [/3dC,™sin(t^Mt-</))|*n+£>(«'„|+li.c.] , (12) 



where the coefficients C™ = (5'„+f |z|'I'„)/d represent 
the real-valued overlap integrals between resonantly- 
coupled WS states. Similarly to the approximation made 
in the AM case, in Eq. ([T2|) we disregarded all off- 
resonance time-dependent coupling terms between dif- 
ferent WS states; in addition, it is worth underlin- 
ing that the diagonal time-dependent terms are here 
kept, in contrast to the AM case, as they exhibit a 
spatial dependence on the lattice index n. Apply- 
ing the transformation into the rotating frame U^f = 
eM-iEn=-^nA{t) l^-n^^ with A(t) = iujM/i)t ■ 
f3d/{hojM) cos{ujMt — 
the time- independent form 



the Hamiltonian takes 



n— — 00 



n— |*„)(*«l + 



+ fe'(-/2-«^|v&„+,)(vI/„|+h.c 
where the tunneling rates J7?^ are defined as 



PM 



iuJB "^^^huM J 



(13) 



(14) 



with Ji the first-order Bessel function of the first kind, 
and ££ the £-th Fourier harmonic of the energy band E{k) 



d 

2^ 



dfc e 



-itkd 



E{k). 



(15) 



The quantity 2 is well known for being the tunneling 
rate between states lying (. sites apart in a pure lattice 
potential, i.e. in the absence of any modulation or ex- 
ternal force. A detailed derivation of the Hamiltonian in 
Eq. (jl3p is given in Appendix |^ 

It is far evident the similarity between pUj) and ([T^]) . 
the only and significative difference being given by the 
actual expression of the tunneling rates ([TT|) and ([T4|. 
In AM case the tunneling rate is linear in the driving 
amplitude a, while in PM case this occurs only for j3d <C 
huiB, when Ji{(3d£/{hujM)) ~ I3d£/ {2hujM)- Besides, in 
the strong-driving regime defined by large values of its 
argument, Ji vanishes at specific values of the argument 
and the tunneling rate is dramatically suppressed leading 
to dynamic localization [H, [H, . 

In both cases, it is demonstrated that the stationary 
states of ([TU| and ([T^ at resonance (i.e, 6 = 0) are delo- 
calized Bloch states with their energies forming a band 



-Emod(fc) = Jesm{ikd- 



(16) 



where Je is given, respectively, by or jf"^ [H, [2] 



1. In real space 

We investigate the time evolution in real space of an 
initially localized state under the action of AM and PM 
driving. The transport dynamics is governed by the 
Hamiltonians in Eqs. ((TU)) and ([T5|) . respectively, which 
are similar in their form. 

For S ^ 0, the translation invariancc is broken and the 
Hamiltonians exhibit a strict analogy with the Hamilto- 
nian of a static lattice in the presence of an effective ho- 
mogenous force of magnitude Fs defined hy 6 = Fgid/h. 
This suggests the physical picture in which the nearly- 
resonant driving reduces much the "tilt" of the lattice 
potential stemming from the force F, but not completely: 
a small detuning S effectively corresponds to a small ex- 
ternal force Fs. One thus expects that the transport dy- 
namics is characterized by Bloch oscillations with time 
period 27r/(5. Indeed, this was first observed through 
macroscopic oscillations of atomic wave packets [25| , and 
later studied with non-interacting Bose-Einstein conden- 
sates (BECs) [4^; it also led to a generalization of Bloch's 



acceleration theorem for periodically driven systems [41 
In the other case with 6 = 0, the system is invariant under 
discrete translations, and coherent delocalization occurs, 
with the atomic wave packets spreading ballistically in 



time according to the dispersion relation in Eq. [24 
The eigenstates of the Hamiltonians are the Wannier- 
Stark states 



|$„(5)) = J_™(^A ) 



rn— — 00 



(17) 



corresponding to a force -F^, with J,„ the m-th order 
Bessel function of the first kind [s^l • These states are cen- 
tered on the n-th lattice site and extend over a number 
of sites of the order of the argument of the Bessel func- 
tions, Je/iFsd). The set of energies associated to these 
states forms a Wannier-Stark ladder with En{S) = nh5/£. 
Hence we obtain the time evolution operator in the ro- 
tating frame 

Urit) = cxp I- if^ n5t/e\<!>n{6)){^n{5)\] . (18) 



By applying it to the initial state \'^n) we get the time- 
evolved state 



|^'„(t)) = e-*"-^*/^ e-™"^e-™''*/2 : 



f smiSt/2) \J,\t\ 

at time t, with the intermediate steps reported in Ap- 
pendix |B] The superposition principle allows the evolu- 
tion of any initial distribution to be calculated, as well. 
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One can use Eq. p9)) to calculate the spatial 
broadening of a matter-wave packet initially localized 
at the n-th site as a function of time, cr^^ (t) = 
(^'„(t)|(z — n(i)^|^„(<)). Considering the expression in 
Eq. ([T9| , one expects the quantity (t) /d to be of or- 
der of the argument of the Bessel functions. This is, in 
fact, confirmed by carrying out the detailed calculation, 
which produces at large t the asymptotic ballistic expan- 
sion 

cr*„ (t) ~ v^t^ smciS/Tpf , (20) 

with Vi = hr^ ld\Ji\l \/2 the speed of broadening at reso- 
nance, Fi? = 27r X (Trf)"""^ the Fourier-limited linewidth of 
the resonance spectrum, and the sine function defined as 
sinc(x) = sin(x)/a;. Appendix [Cl reports a derivation of 
this expression for the resonant case ((5 = 0). We omitted 
showing in Eq. (|20p the terms corresponding to the size 
of the initial state |4'n) and a quivering motion exhibit- 
ing fast, tiny oscillations on a timescale of the order of 
2-K /ujm] in fact, these terms are of marginal significance, 
as they depend specifically on the given initial state. The 
complete expression is however provided in Appendix [U] 
Thus, from the expression ([^0)1 we see that resonant 
driving has the effect of coherently delocalizing atomic 
wave packets over sites which arc distant id apart; con- 
versely, nearly-resonant driving yields spatial oscillations 
with period 2n/5 and maximal extent 2vi/5. The pre- 
diction in Eq. (|20p for (T*„(i) allows the measurement 
of Bloch frequency by directly determining the resonant 
frequency (.ujb through in-situ observations of the spatial 
distribution of the atoms in the lattice. In reality, one has 
to consider that the atomic cloud is characterized by an 
initial finite size, ao, which produces a spatial distribu- 
tion that is the convolution of the two distributions: the 
variance, corrected for the atomic cloud's size, becomes 



a{t)^^ai + v^,t^.inc(^^^^i^y. (21) 

This expression depends on four free-parameters, 
which can be determined by fitting the experimental 
data: the starting width of the atomic cloud ao, the max- 
imal extent at resonance vgt, the resonance linewidth F, 
and the Bloch frequency lob- The experimental valida- 
tion of this formula will be analyzed in Sec. IIVI 

2. In momentum space 

Studying nearly-resonant driving in momentum space 
brings an alternative viewpoint, which puts in the fore- 
ground mainly the aspects related to quantum transport. 
One in fact expects that at resonance the eigenstates of 
Eq. (|18p are delocalized Bloch states \k) with wavevector 
fc, and the transport dynamics is determined by the en- 
ergy dispersion relation E^odik) introduced in Eq. p6l) . 



Conversely, when off resonance the eigenstates are the 
Wannier-Stark states |$„((5)) given in Eq. ([T7|. and one 
expects that the transport-dynamics is characterized by 
Bloch-like oscillations: the quasi-momentum increases 
linearly in time as hdk/dt = —F^t, with the understand- 
ing that k is taken modulus the Brillouin zone. 

This can be directly seen by projecting the time- 
evolved state ^„(i)) on quasi-momentum states, 

(fc|M/„(t)) = (fc|M/„) e-™^*/^x 

X exp dt' £;„,od(fc + Fst'/h)^ , (22) 

where we made use of Bloch's theorem which provides 
translational symmetry to the WS states 

(fc|vI/„+„) =e-""^'^(fc|vI/„) (23) 

and the Jacobi- Anger expansion to derive it. 

This result allows us to interpret an optical lattice that 
is dynamically driven by AM or PM modulation as a 
tilted stationary lattice whose parameters can be tailored 
at will. The tilt is controlled though the detuning 5. The 
periodicity of En\od{k) as a function of k shows that the 
size of the Brillouin zone is shrunk by the order of the 
driving harmonic t. This is rigorously proven in the next 
section. The overall phase of the effective energy band 
(j) can be tuned through the phase of the time-dependent 
external driving. The width of the energy band Je, can 
be manipulated as well by tuning the driving strength 
a and /?. In the PM case the width depends also on 
the ratio /3£d/{hu!M), due to the mechanism of dynamic 
localization for strong driving. 

Besides, the result implies that the momentum 
spectrum of an initially localized WS state evolves in 
time in the presence of external driving by only acquiring 
a phase, which depends on the time and on the wave 
vector k. Then, the squared modulus of the momentum 
spectrum remains constant 

\{k\^,M\'^\{k\^o)\\ (24) 

implying that only the probability amplitudes and not 
the probability distribution in the momentum space 
evolve in time. 

However, dynamical driving can be used to extend spa- 
tial coherences of wave packets, i.e. delocalizing them 
over several lattice sites, leading to transient narrow dis- 
tributions in time-of-flight (TOF) experiments. Tran- 
sient modification of TOF density distributions is at the 
basis of DEBO mechanism, and it is discussed in Sec. 
Ill El in static optical lattices, transient interference phe- 
nomena have been recently investigated using BECs [42| . 

D. Transport resonance spectrum 

We want to study the problem of quantum transport 
in driven optical lattices from the point of view of space 
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and time symmetries. In particular, we are interested 
in determining the exact positions of resonances without 
making any approximation of the original Hamiltonian in 
Eq. ©. This is of special interest for the application of 
AM or PM driven quantum transport to determine the 
Bloch frequency cjs in Eq. (j4|) , thus performing a precise 
and accurate measurement of the local force F. 

It is convenient to define 'Ho as the Hamiltonian in the 
absence of modulation, i.e., for a = 0, /? = 0, and Hi as 
the driving Hamiltonian Hi = % — Ho- Regarding space 
symmetries, we have that Dm'Ho'Djj^ = T-Lq + mhujB and 
Vm'Hi'D^n = Hi + mf3df{t) with Dm the unitary oper- 
ator shifting the space coordinate by an integer number 
TO of lattice sites. Regarding time symmetries, we have 
that Ho docs not depend on time, while Hi is invariant 
under time shifts by an integer number of tm = 2tt/ujm- 

The invariance under discrete-time shifts, which occurs 
for periodic drivings, allows the use of the Floquet theory 
[26j in order to express the time evolution operator lATit) 
at the time t = k tm as 

UT{tiTM)^UT{TMT (25) 

with K an integer value. The Floquet quasi-cnergies 
Ei are determined, up to integer multiples of 2t:ti^[ /ft, 
by means of the cigenstatcs cxp{'~i EiTM / fi) of IAt{tm): 
transport resonances manifest themselves when the 
corresponding eigenstates of Ut{tm) are delocahzed 
states, leading to coherent delocalization of any initially- 
localized wave packet. In addition, as a result of Bloch's 
theorem, we have that the eigenstates are indeed delocal- 
ized when 'DmUT{TM)T^m = ^t(ta/), i.e., when Ut{tm) 
is invariant under translations of an integer number of 
lattice sites. 

By making use of the interaction picture, the evolution 
operator can be written as 

Ut{tm) = cxp{—iHoTM/h) Texp ^— z J dtHi/hj , 

(26) 

where Texp represents the time-ordered exponen- 
tial and Hi — cxp{iHoTM /fi)HiCxp{—iHoTM/fi) 
is the transformed Hamiltonian. Owing to the 
space symmetry mentioned above, it follows that the 
first term transforms as I?„i exp(— z'HoTA//ft)2?J„ = 
exp{—iHoTM / fi) cxp(— z?7iwbtj\/), while the second one 
remains invariant 

Texp (-i J dt Hi/hj X>t„ = 

= cxp (-ij^ dt ml3df{t)/hj Texp (-i "dtHi/hj , 

(27) 

where the identity J^^'dtf{t) = is enforced by the def- 
inition of Tm- It results finally that 

I?„Z^T(TM)2?r„ ^ UtItm) cxp{-imujBTM), (28) 



implying that [Z/r(rj\/), Dm] = for ujm = ujb/q with 
q an integer number. More generally, if we consider the 
symmetry group characterized by discrete translations of 
£m lattice sites, with £ and to integer numbers, it fol- 
lows that [Ut{tm), Dim] = for lom = i'-^s/q, proving 
that transport resonances, i.e. coherent delocalization of 
atomic wave packets, occur at integer (and sub integer) 
harmonics of the Bloch frequency ub- 

In addition, this result proves that, when modulat- 
ing at integer harmonics the translation symmetry 
group is the one determined by a unit step of £ lattice 
sites, indicating that the first Brillouin zone which de- 
scribes the eigenstates oiUriTM) extends over the inter- 
val [— 7r/(£d), 7r/(£(i)], thus showing indeed a shrinking 
by an integer factor £ of the Brillouin zone of the original 
lattice. 

These findings can be straightforwardly generalized to 
an arbitrary periodic driving f{t) with period lixji^M- 
It results that for AM driving the location of transport 
resonances remains unchanged, while for PM driving one 
should consider that a non-zero average value of the driv- 
ing force F' = J^" dt P f{t)/TM produces an effective 
Bloch frequency shift F'd/h with respect to ujb defined in 
Eq. |4l This generalization to arbitrary periodic drivings 
is of clear interest for measuring ujb with high precision 
and accuracy: measurements obtained by AM driving 
are insensitive to the presence of spurious higher and sub 
harmonics of w^/ hi the modulation signal f{t). 

E. Delocalization-enhanced Bloch oscillations 

In Sec. Ill C II we showed how resonant driving of the 
lattice potential establishes coherences among WS states 
as a result of the broadening of atomic wave packets dur- 
ing the modulation time. In the case of a finite time of 
free evolution after releasing from the optical lattice, this 
implies interference of matter waves initially located at 
different lattice sites. This can be put into evidence by 
means of the following three-steps experiment: we apply 
the (AM or PM) modulation f{t) = sm{ujBt) for a certain 
time tjnod, then hold the wave packets in the static lattice 
for thoid, and finally release the wave packets in time of 
flight ttof- For an optimal choice of tmod and ttof it is 
possible to enhance the interference effect responsible for 
BOs and its visibility. 

To understand the DEBO mechanism, we thus start 
with an initially localized state jV'o) centered around a 
site n = 0. The results could be easily generalized to any 
initial site n by a simple translation. The initial state 
at t = can be written as in If the thermal energy 
/csT is larger than the recoil energy Efi, G{k) is smooth 
and flat over the whole Brillouin zone. In the first step, 
after a modulation at the ^-th harmonic for a time tmod, 
the states \ipB{k)) thus evolve by means of the operator 
U{t) = cxp{-iEmod{k)t/h): 

\Mk,Unod)) = e-^^^'"(^'^-'')|V'B(fc)) (29) 



where the quantity A = Jn tmod/f^ multiphed by i repre- 
sents the spatial extent of the broadened wave function 
expressed in units of the lattice constant d. 

During the hold time thoid in the static lattice, the 
external force shifts k according to A; — ^ fc + coBthoid/d 
and the wave function changes into 



li'itmod + thold)) = 

= / dkG{k)e-'^''''^^"''''>\Mk + iOBtkoid/d)) (30) 



The result ([50)1 can be alternatively stated as: 




\i^{t,nod + thold)) 



dkG{k - iOBthoid/d)x 



where the argument of the G function has to be taken 
modulus the Brillouin zone, and tb is the Bloch period. 
Switching off the lattice maps the Bloch states into free 
particle states. This procedure can occur suddenly or 
adiabatically. While in the first case complex interfer- 
ence effects involving higher energy bands are expected 
to appear, we want here to catch the basic concept in the 
simplest case of adiabatic release. After a time of flight 
ttof , the probability amplitude to find the atom along the 
optical lattice direction is: 



{zlll'itmod + thold + ttof)) = / dkG{k-UBthold/d)x 

Jbz 



The integral in ([5^ determines the shape and posi- 
tion of the wave packet. In particular we are interested 
in the interference peak which is produced as a caus- 
tic effect because of dispersionless evolution at the in- 
terference points. Since we supposed a flat distribution 
G(fc), the interference peak is explained in terms of a 
vanishing second derivative of the phase in the exponent 
g{k) = A siniikd - 2TTethoid/TB) + hkHtof/i2M). We 
also require a vanishing third derivative to make the re- 
gion with g"{k) = as large as possible. The two condi- 
tions are satisfied when ttof = Ttrans where 



Ttrans = MA^d'^/h 



(33) 



is the transient time after which the visibility of the in- 
terference peak is maximal. Fig. [T] clearly displays the 
emergence of a dispersionless peak in the wave packet as 
the time of flight approaches Ttrans- 

The transient time Ttrans (|33p has a clear physical 
explanation. Since fc^T > En, the wave packet is es- 
sentially localized in a single lattice site, from where it 
broadens in time of flight with a speed h/{M£d). Only 



FIG. 1. Enhancement of the interference peak during time- 
of-flight detection of Bloch oscillations for a thermal wave 
packet {kBT/Eii ~ 3) by means of a resonant modulation 
burst [l = 1). The solid (blue) line represents the single 
atom wave function after the three steps BO measurement 
with a broadening factor yl = 30 lattice sites and ttof = 
MAd'^/{2h) which exhibits a narrow feature at the center 
of the wave packet with respect to the wave packet without 
initial broadening (blue dashed line). The third (black dash- 
dotted) line shows the wave function at ttof = M Ad^ / {h) = 
Ttrans which shows the narrowest peak. 



after a time Ttrans the time-of-flight expansion is equal 
to the extent of the broadened wave packets after modu- 
lation which is id A. Picturing the wave packet as a su- 
perposition of different waves initially located in different 
lattice sites, the maximal interference occurs when all the 
broadening waves in time of flight meet each other, which 
occurs at the transient time Ttrans- The DEBO technique 
can be applied for measurements with superior sensitiv- 
ity to the effect of external forces. The details of the 
experimental results are presented in Sec lIVI 



F. Decoherence in momentum space and resonant 
tunneling 

Since both BO and WS resonant tunneling are co- 
herent quantum phenomena, it is important to investi- 
gate the effects of decoherence. Basic decoherence mech- 
anisms are spontaneous emission, atom-atom interac- 
tions [1^, lattice random recoils [3] and environmen- 
tal scattering (i.e., off- resonance scattering by incoherent 
light or background gas collisions). 

While both spontaneous emission and environmental 
scattering do not affect the dynamics of an ultra-cold 
gases in a high vacuum system and in their ground state, 
it is possible to evaluate the effect of atom or photon 
scattering. They can in fact have different strengths, 
leading to a collapse of the atomic wave function in a 
single site, thus leading to a decrease of the BO visibility 
in momentum space or sequential tunneling in position 
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space, even in the presence of a modulation. 

The effect of a decoherence process in momentum space 
can be understood in terms of superposition of WS states. 
Starting from ([8]) (with a = /3 = 0) and exploiting the 
translational symmetry of WS states (P^ . we can rewrite 
a general expression for the atomic wavefunction in the 
momentum space (45j 

oo 
n— — OO 

= {k\^o{k)) C{k + Ft/{dh)) . (34) 

Thus the wave function is the result of the product 
of a single WS momentum distribution and an envelope 
function C(fc), which is a comb-like distribution periodic 
both in time and in momentum space. The width of 
the momentum peaks AA: is inversely proportional to the 
number TV of the superimposed WS states, as we already 
showed in the previous section. Any scattering event like 
photon recoil or atom-atom collisions yields a decay of the 
WS states superposition, while conversely a diffusion-like 
behavior of the wave packet dispersion. Thus we expect 
that the momentum peak width evolves with time as 




FIG. 2. Experimental setup. The output radiation from a sin- 
gle mode 532 nm (Verdi) is split and then focused on the ultra- 
cold **Sr cloud, while a pick-off is sent through 200 m fiber 
to the comb lab for frequency measurements, fbv- beat-note 
frequency between Verdi and comb; /q: the carrier-envelope 
offset frequency, frep' repetition rate frequency; SMF: single 
mode fiber; AOM: acousto-optic modulator. 



Afc2(i)~ Afc(0)2(l+7t) 



(35) 



where 7 is the rate of the scattering process. For an exact 
treatment of the dynamics of the wave-packet, and thus 
of the peak width Afc, it is necessary to solve the stochas- 
tic Liouville equation (46j . where a localization operator 
is introduced and a master equation for the evolution of 
the density matrix p is provided. The prediction for the 
momentum peak broadening in psp is equivalent to say 
that a random Gaussian process heats the atomic gas un- 
til a uniform distribution of the atoms over the Brillouin 
zone is obtained. This random walk in the momentum 
space is proportional to the time t and the density of 
scatterers (the atoms or the lattice photons). 

In the case of BO experiment, wc must also con- 
sider the presence of the excited bands which introduce 
additional source of losses. In fact for shallow traps 
(C^o ^ Eft) collisions can result in the transition of the 
atom from the lowest band into the upper bands which 
are only weakly bounded. This means that only collisions 
with small exchanges of momentum can effectively con- 
tribute to the broadening of Afc and hence ^exp < 7 |44 1 . 



The lattice-recoil effect plays a smaller role since atom- 
lattice scattering energy exchange excites the atoms into 
the second band (the energy gap between the first and 
the second band is typically Eq ~ En) and thus does 
not contribute to the BO signal. 

In the case of resonant driving applied to the system, 
the localization process inhibits the coherent coupling 
among adjacent WS states which is responsible of the 
resonant tunneling of the initial atom distribution. The 
decay of the off-diagonal elements of the density matrix 



causes a sequential tunneling which asymptotically gives, 
as already stated, cr^(t) Dt. The sequential tunnel- 
ing processes have the effect of broadening the resonance 
spectrum. In the case of simple localization Lindblad op- 
erators Ln = ^\n){n\^ the numerical integration of the 
master equation leads to a Lorentzian-shaped resonant 
tunneling spectrum with 7 approximatively the linewidth 
of the spectrum [l^] . The presence of a decoherence pro- 
cess will result into a broadening of the Fourier-limited 
linewidth nearly by an additive quantity 7, which is the 
inverse of the sensor coherence time. 

Having presented decoherence from both points of view 
of momentum and real space, we will discuss quantita- 
tively the impact of decoherence on the final sensitivity 
in Sec. El 



III. EXPERIMENTAL SETUP AND 
PROCEDURES 

The experimental setup is based on cooled and trapped 
^*Sr atoms The details of the atomic sample prepa- 
ration was already described elsewhere [i^. Here we 
limit our description to those basic elements which have 
a relevance on explaining the results in terms of accuracy 
and sensitivity of our atomic force sensor. A schematic 
picture of the experimental setup is shown in Figj^l 

Atoms from a thermal beam are slowed by a Zee- 
man slower and trapped in a "blue" magneto-optical trap 
(MOT) operating on the Pi resonance transition 

at 461 nm. The temperature is further reduced by a 
second cooling stage in a "red" MOT operating on the 
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^5*0 Pi intercombination transition at 689 nm. This 
produces about 10^ atoms at a temperature of 1 ^K. 
Since the force of gravity is comparable to the force ex- 
erted by the red MOT on the atoms, the cloud of trapped 
atoms assumes a disk-like shape with a vertical size of 27 
/im and a radial size of 180 /im. The atoms are adi- 
abatically loaded in an optical lattice in 200 /is. The 
lattice potential is generated by a single-mode frequency- 
doubled Nd:YV04 commercial laser (Verdi- V5 Coherent, 
Xl = 532 nm) delivering about 1 W on the atoms. The 
optical lattice is then produced either by retro-reflecting 
the input beam or by means of a second fiber which de- 
livers half of the available power from the bottom side of 
the vacuum chamber [i^ . A telescope can set the beam 
waist between 300 and 560 ^m at will. 

Typical experimental parameters yield a lattice popu- 
lation of Nat ~ 5 • 10*, for an atomic density n = 6 ■ 10^^ 
atoms/cm"^ and lattice depth Uq = 3Er. If wc con- 
sider the decoherence due to elastic collisions driven 
by a scattering length ass = — 1.6ao [EO], we expect 
7co« < 2.2 • 10""^ s~^. A high vacuum system allows 
to exploit the long coherence time allowed by strontium 
atoms. It consists of a 50 l/s ion pump and a titanium- 
sublimation pump yielding a background pressure level 
smaller than 10"^'^ Torr. The measured lifetime in the 
vertical lattice is about 20 s. 

We used earth gravity as a test force for the atomic 
sensor. This is possible by placing the optical axis of 
our lattice along the direction of g. The beam is aligned 
along the gravity axis as follows: the superposition of one 
beam onto the other is obtained by maximizing the trans- 
mission through the fiber of the opposite beam. Then 
the lattice verticality is made sure by interposing a large 
glass container with water over the vacuum cell where 
atoms arc trapped. Part of the upwards-directed beam 
is reflected by this surface with an angle 9. Finally we 
measure far apart (£> = 3 rn) from the cell the distance 
s ~ 6D between the laser spot formed by this reflection 
and the second downwards-directed laser beam passing 
through the water container. This distance s is then re- 
duced to zero through the fine adjustment of the optical 
table tilt. This experimental procedure then prevents 
any possible systematic shift due to both relative and 
vertical alignment of the two counter-propagating laser 
beams. The residual systematic uncertainty is described 
later. 

Absorption imaging of the atomic cloud position and 
distribution is performed either in situ or by means of 
time-of-flight (TOF) technique by a CCD camera with a 
spatial resolution of 4.6 /im. 

Different procedures can be applied to measure g. In 
the case of PM resonant tunneling, a piezoelectric ac- 
tuator was attached to the retro-reflecting mirror. The 
AM technique, instead, employs one double-pass acousto- 
optic modulator (AOM) which is servo controlled to si- 
multaneously stabilize the lattice depth and provide the 
modulation frequency ujm to the lattice. This modula- 
tion is either continuous or with a sequence of modula- 
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FIG. 3. Noise spectra of the relative intensity of the lattice 
laser at 532 nm which is used to produce the optical lattice po- 
tential. An intensity stabilization system strongly suppresses 
the intensity noise of the laser. The comparison of the two 
spectra shows that the servo lock attenuates the bump of noise 
originating from the acoustic vibrations by about 30-^40 dB. 



tion bursts depending on the kind of experiment. The 
former operational method is used to set resonant tun- 
neling, the latter for DEBO or Loschmidt echo (which 
is not discussed in this work, but detailed in Ref. [23| ) 
experiments. Both the modulation frequency and the 
experimental timing are referenced to a GPS-stabilized 
signal, as shown in Figl2l 

A servo control system is employed to stabilize the laser 
intensity, and thus avoiding either slow drifts of the lat- 
tice depth Uo and suppressing acoustic vibration-induced 
amplitude noise. The laser intensity is controlled by 
changing the radio frequency (RF) signal which drives an 
AOM. The error signal is then obtained by comparing the 
intensity measured from a pick-off after the optical fiber 
and a stable voltage reference. Hence, any misalignment 
of the fiber injection does not affect the lattice potential 
depth. Instead fluctuations of the polarization after the 
optical flber can also affect the potential depth, mostly in 
a conflguration of two counter-propagating laser beams 
with respect to retro-reflected single beam. The stabi- 
lized intensity noise spectra in Fig. [3] shows the reduction 
of the free-running amplitude noise in the acoustic range 
by about two orders of magnitude, with a servo band- 
width of about 20 kHz. The noise frequencies which are 
most harmful to the atoms are, in fact, in the range from 
Hz to several kHz, which effectively coincide with the typ- 
ical acoustic frequencies. Excitations at these frequencies 
induce heating of the atoms by exciting the vibrational 
frequencies of the optical lattice [5l|. As shown in Fig. |31 
we achieve a stabilization of the intensity offset at the 
relative level of 10^®. 

The servo loop controlling the lattice laser beam in- 
tensity also allows to insert the AM driving signal which 
shares the same frequency reference as the other relevant 
frequencies. From the recorded FFT spectra of the mod- 
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ulated laser intensity signal the higher harmonics con- 
tents have been reduced of more than 60 dB, i.e. below 
the noise floor level. This reduces the projected Bloch 
frequency uncertainty at 20 ppb level (see Sec. IV Bp . It 
is evident that the possibility to have an in-loop driving 
modulation in the AM resonant tunneling and a direct 
measurement of the driving into the lattice make the AM 
technique more suitable for high precision measurement 
with respect to the PM technique employing a mechani- 
cal transducer. 

The PM technique is performed by applying the driv- 
ing signal to the PZT attached to the retro-reflecting 
mirror. The response of this actuator was characterized 
by an heterodyne measurement with a Michelson inter- 
ferometer, which shows a linear dynamical response of 
the electro-mechanical system with a cut-off frequency of 
4 kHz [13]. The voltage applied to the PZT allows to 
modulate the position the lattice potential by up to 10 
lattice sites peak to peak, so that zq < 2.5Xl- According 
to (|14p . this means that the non-linear regime showing 
dynamic localization can occur only at modulation fre- 
quencies higher than the 6-th harmonic. 
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FIG. 4. Example of frequency counting of Verdi laser for 
calibration of the ID **Sr optical lattice frequency, reporting 
for comparison the typical duration of gravity measurements 
for both techniques. The beat-note signal between lattice 
laser and the optical comb is used to stabilize the repetition 
rate of the femtosecond laser. The intervals where no data is 
reported are caused by the optical frequency comb running 
out of lock from the lattice signal. 



A. Calibration of the lattice frequency 

A precise knowledge of the lattice laser frequency value 
and its instability is needed for high precision force mea- 
surements. We employed a Ti:Sapphire femtosecond fre- 
quency comb to obtain the absolute frequency measure- 
ments of the lattice laser and to characterize its frequency 
stability [52| . As shown in the schematic diagram of the 
experimental setup in Figj^J the output radiation from 
the 200-meter-long fiber is superimposed with radiation 
from the femtosecond comb (filtered at A = 532 ± 5 nm) 
and the beat signal with the corresponding tooth is ob- 
served on a fast photodetector. 

To perform precise frequency measurement, we lock 
the repetition rate of the comb by stabilizing this beat- 
note to a RF-frcquency synthesizer and we count the 
carrier-envelope offset frequency of the optical frequency 
comb by self-referencing techniques. The absolute value 
of the lattice laser frequency is given by equation / = 
n X frep ± /o i where the correct tooth number n 
is determined by a wavemeter calibrated with a stabi- 
lized laser on the 689 nm ^Sq-^Pi ^*Sr transition, that 
allows an evaluation of the lattice laser frequency with 
an uncertainty of the order of 100 MHz. 

Fig 13] displays the result of an absolute frequency mea- 
surement that lasts for about 10^ s. It is interesting to 
compare dynamics of the lattice laser frequency with typ- 
ical durations of Bloch frequency measurement both with 
BO technique and with AM resonant tunneling. From 
this measurement we can clearly observe periodical fluc- 
tuations at two different regimes: the first low frequency 
fluctuation with typical period of about 10^ s and am- 
plitude of 130-140 VIHz, and a second regime is a fast 
oscillation with a shorter period of about 33 s and a typ- 



ical amplitude of 30-70 MHz (see inset in FigH]). 

Long term deviation of lattice laser frequency greatly 
depends on many factors, mainly related to temperature 
and laser power setting. To increase the long term sta- 
bility we left the laser continuously operating for about 
two days at the same power setting. Moreover, a thermal 
insulation foam is applied to the laser head to increase 
the passive temperature stability. However, this model of 
commercial laser can show both quite different absolute 
frequency (offset by several GHz) and large frequency 
drift rate at start up (up to 100 MHz in one minute). 

Hence, long-lasting measurements of forces by means 
of the Bloch frequency are affected by lattice frequency 
uncertainty at the 100 ppb level. Looking at the length of 
the two force measurements, BO measurements usually 
requires 5.6 • 10^ s (1.5 hours), while AM resonant tun- 
neling lasts three times less. Then the sensitivity of BO 
measurements is three times higher, increasing the er- 
ror for each Bloch frequency determination, though the 
systematic effect can be equivalent in both techniques. 

IV. MEASURING EXTERNAL FORCES: 
BLOCH OSCILLATIONS VS. DRIVEN 
RESONANT TUNNELING 

A. Effects at resonant driving I: Wannier-Stark 
localization vs. coherent delocalization 

As discussed in Sec. Ill CTl atomic wave packets loaded 
into a driven vertical optical lattice potential exhibit co- 
herent delocalization arising from intraband transitions 
among WS levels. WS intraband transitions arc ob- 
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served by monitoring the in-situ wave packet extent. Co- 
herent delocahzation sets in for modulation frequencies 
lom = lob, or multiple integers € of a;^, suggesting that 
tunneling occurs not only between nearest neighboring 
sites (i.e., when £ = 1) but also between sites that are 
I lattice periods apart. We assessed this phenomenon by 
observing an increase of the atomic distribution width 
following Eq. (PT|) . 

In Fig. [5] the effect of AM driving is shown: after res- 
onant amplitude modulation at frequency ojm = 6 x wb, 
the recorded density profile of the atomic cloud appears 
clearly stretched along the lattice direction, while the 
density height decreases. 

The measurement of the spatial rms extent along 
the lattice for several modulation frequencies allows to 
determine the Bloch frequency with high sensitivity. 
Long modulation times (tmod = 10 s) allow to per- 
form intraband WS states spectroscopy with very narrow 
linewidths, down to Fourier-limited linewidth Tp = 31 
mHz. An AM resonance spectrum for the 6-th harmonic 
is shown in the inset of FiglSj Here the observed F corre- 
sponds within the experimental error to the Fourier lim- 
ited linewidth Tp- Long modulation times also increase 
the transport resonance quality factor as the maximal 
extent at resonance is vetmod- However, as observed in 
Fig. \E\ this implies a reduction of the observed optical 
density and thus a reduced signal-to-noise ratio of the 
in-situ absorption imaging of the atomic cloud. The nec- 
essary trade-off between an enhanced quality factor and 
a reduced signal-to-noisc ratio is found by an optimal 
choice of the modulation depth a. A detailed analysis of 
the effects which may limit the WS resonant tunneling 
sensitivity is presented in Sec. IV Bl 

Comparable results were previously achieved with the 
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FIG. 5. Effect of the resonant AM technique on the atomic 
cloud with a modulation time of 10 s. The integrated density 
distributions after AM for different detunings 5 from 6xajs 
clearly appear stretched at the resonance. In the inset, we 
plot the resonance spectrum derived from the rms size along 
the lattice for AM resonant tunneling at the 6-th harmonic. 



PM resonant tunneliiig. Coherent delocahzation was ob- 
served up to £ = 4 [2J] and Fourier-limited linewidths 
up to 15 s. In this case the existence of higher order 
harmonics in the lattice potential is unavoidable. The 
measurements show in the PM case a higher fluctuation 
of the broadening velocity ve which can be attributed to 
the non linear contribution of the higher order harmonics 
as function of the modulation strength /3. For this rea- 
son and according to possible frequency shifts presented 
in Sec. Ill D( from now on we will consider only AM reso- 
nant tunneling technique as an accurate sensor for force 
measurements . 



B. Effects at resonant driving II: 
DEBO experimental results 

We experimentally compared the DEBO technique 
with the usual BO measurement for high sensitivity grav- 
ity determination. Enhancement of the interference peak 
during timc-of-flight detection of BOs was previously 
used for coherent control of the spatial extent of an 
atomic wavefunction [2^ for short times. We applied 
more recently this technique for high sensitivity force 
measurements [H- 

The phase sensitivity to BOs is greatly enhanced by 
means of DEBO described in Sec lIIEl In order to demon- 
strate this, we employ the following sequence: we adi- 
abatically load the atoms form the red MOT into the 
shallow lattice {Uq = 3Eji) so that after few milliseconds 
the in-lattice temperature is reduced to about 0.6 fiK. 
We then apply an AM burst at the Bloch frequency last- 
ing 120 cycles with a modulation depth a = 0.2. After 
different time intervals t of free evolution in the lattice, 
we perform absorption imaging of the atomic cloud by 
means of TOF technique and the resulting momentum 
distribution is analyzed. 

The result of the DEBO technique is shown in FigEl 
we are able to follow BOs with high sensitivity due to 
the formation of narrow interference momentum peak, 
as shown in the left inset in Fig. |6l The resulting in- 
terference pattern shows a narrow Gaussian distribution 
due to the convolution of the wavefunction shown in Fig[T] 
with the initial atomic spatial distribution. The effective 
momentum peak width is Ap sa 0.18 fifc/, for a ttof = 14 
ms, which corresponds to Ttrans in p3|) . This result is 
nearly a factor of two narrower than the central momen- 
tum peak in non-interacting BECs [s^ [sJl for the same 
time ttof- 

The comparison with BO without initial AM burst is 
also shown in Fig|6] (blue circles). In this case, the in- 
terference pattern is described by two overlapping Gaus- 
sian distributions overfilling the Brillouin zone [ll| . The 
recorded single point error ap can be as large as 0.2 ^/cl, 
for an overall Bloch frequency relative error Auje/i-ob = 
3 • 10^^. The DEBO interference peak yields a typical 
error of each single point of the momentum evolution 
(cTp) = 1.6 • 10~^ftfci. This result is due to the greater 
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FIG. 6. Comparison between Bloch oscillation phase evolu- 
tion with initial AM burst (black diamonds) and without it 
(blue circles). In the insets, 2D TOF atomic distributions are 
displayed for the two cases. The red line is obtained from the 
best fit with a sawtooth function. 
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FIG. 7. Summary of gravity measurements by means of 
DEBO (black squares) and AM resonant tunneling (red cir- 
cles). The dashed line represents the precise value of the 
gravity g measured by the FG5 [H which we use to compare 
the two force measurement techniques. 



visibility obtained by the AM burst tcclinique: the peak 
visibility, defined as in Ref. [42|, can be as large as 60%, 
while for classical BOs it is smaller than 10%. Implica- 
tions on the precision in determining ujb are discussed in 
Sec. TO 



FORCE SENSITIVITY AND SYSTEMATICS 
OF THE ATOMIC SENSOR 



We analyze the sensitivity and the final precision on 
the determination of lub and its application on the de- 
termination of the gravity acceleration g. We proceed 
by summarizing the systematic effects which affect the 
system for both techniques. 



A. Phase sensitivity in Bloch oscillations 

Starting from an almost dispersionless momentum dis- 
tribution due to the AM resonant burst, the increase of 
visibility yields the observation of BOs up to 17 s, mainly 
limited by the lifetime of the optical trap. Each gravity 
determination lasts more than 2 hours. Repeated mea- 
surements with the same experimental conditions as in 
Ref. [l| are shown in FiglT] The precision in the determi- 
nation of the Bloch frequency Alob/i-^b ranges between 
2.2 -10"^ and 1.7- 10~^ with an observation time of about 
13 s. However we observed a relative scatter of the calcu- 
lated gravity values one order of magnitude higher. The 
final result of local gravity is gBO = 9.80488(6) m/s^. By 
looking at the residuals produced by the best fit with a 
sawtooth wave, we observed a reduction of its mean value 
by removing the data set with longer evolution time. This 



is a clear signal that phase sensitivity to BO is washed 
out by long term drift of the apparatus. In particular, 
the main sources of instability which leads to this high 
uncertainty on the determination of gravity are the fre- 
quency of the lattice laser and the MOT center position 
stability, which comes from the long-term instability of 
the second-stage cooling laser. 

It is possible to determine the coherence time of the 
system by looking at the momentum peak broadening. 
According to Eq. ([55]) . if the free evolution time in the 
lattice thoid ^ 7^^ J we can assume a linear peak broad- 
ening Ap{t) — Apo(l + jt/2). In fact, if we consider the 
decoherence due to elastic collisions, the expected coher- 
ence time can be estimated between 440 and 880 s for a 
50% fluctuation range over the number of atom. We ex- 
perimentally observed a slow increase of the interference 
peak width along the 13 s of measurement time. The 
estimated coherence time, defined as the inverse of the 
slope of a linear fit, is equal to 530 ± 150 s, which agrees 
with the expected scattering rate, as resulting from the 
data in Fig. [Sj 

A fundamental source of uncertainty on the phase of 
the BO measurement is the stochastic nature of the os- 
cillation phase due to LZ tunneling during the switch- 
off of the lattice in a time-of-fiight measurement jssj . 
In fact, when the atom cloud is released by the optical 
trap, there is a finite time interval AT^z in which the 
LZ tunneling probability is not negligible. In our exper- 
iment the switch-off time constant is about 5 /is, which 
roughly correspond to a non-zero LZ tunneling proba- 
bility time interval ATlz = 14 ^s and thus a random 
phase 5(j) < ATlz^b affecting the BO of each atom. LZ 
tunneling is then translated into a source of statistical 
uncertainty which can be numerically modeled, for in- 
stance, by means of a Poisson-like Atj^z distribution. For 
our typical values of the lattice depth {Uq < QEb) and 
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atom number {Nat — 5 • 10^), LZ tunneling sensitivity 
corresponds to Soji^z/^b < 10^^. 

B. Sensitivity in Wannier-Stark resonant tunneling 

High sensitivity measurement of resonant tunneling in 
an optical lattice requires coherence between the WS 
states to keep at a very low level any diffusive effect on 
the lattice which can randomly extend the atom cloud 
a thus causing a broadening of the tunneling frequency 
resonance T. 

The keys to enhance the force measurement sensitiv- 
ity are shown in Eq. (|2T|). In the case of Fourier-limited 
resonance (23 |. the lob sensitivity to spatial width fluc- 
tuations Act can be estimated from a power expansion 
of this formula and then by applying a first-order prop- 
agation of uncertainty theory. This results on a Bloch 
frequency uncertainty 

3 

AuJB = 27r X Act . (36) 

nt'^ Vi I 

It is noteworthy that the resonance width Tp = 27t x 
(Trt)~^ does not depend on the index £ which identifies the 
modulation harmonic with respect to the fundamental 
frequency lob- This fact is significant because, as shown 
by expression p6p for the uncertainty on , it allows to 
increase the sensitivity of the measurement by a factor i. 

A typical Fourier-limited AM resonance at the 6-th 
harmonic has been presented in Sec. lIVl in Fig.O In order 
to investigate the sensitivity of the AM resonant tunnel- 
ing measurement, we performed repeated measurements 
of the Bloch frequency with a modulation time t = 10 
s at ^ = 5. Typical uncertainty of the Bloch frequency 



0.22 - 




5 10 15 

Time (s) 

FIG. 8. Broadening of the interference peak in time-of-flight 
measurement as function of the free evolution time inside the 
lattice. From the linear fit, a low density strontium cold gas 
can achieve a coherence time higher than 500 s. 




resonance width difference T - F^ (mHz) 

FIG. 9. Study of the broadening of the AM tunneling reso- 
nance. Considering a modulation time of 10 s, which corre- 
sponds to a relative Fourier linewidth of 31 mHz, it is possible 
to estimate the coherence time of the system by measuring the 
deviation of the experimental linewidth from its fundamental 
limit. 



ranges between 0.1 and 0.3 mHz, which is translated into 
a relative uncertainty between 180 and 590 ppb. 

The statistical uncertainty on the WS tunneling reso- 
nance can be estimated from ()36p by looking at the fiuc- 
tuation of the atomic cloud during the cooling and trap- 
ping procedure. If one assumes an average fiuctuation 
Aa/ao ~ 2%, which is a typical value for our apparatus, 
we get an expected relative error on the Bloch frequency 
of 370 ppb. This effect represents the main source of 
statistical uncertainty for the WS resonant tunneling. It 
can be reduced by starting atomic cloud fluctuations by 
means of a narrow optical tweezer [i^ . 

Any source of broadening of the resonance is trans- 
lated into a smaller quality factor and thus into a reduced 
sensitivity to the external force sensed by the atoms. In 
resonant tunneling experiments, collisions lead to sequen- 
tial tunneling instead of coherent tunneling. Sequential 
tunneling is expected to increase the width of resonance 
spectra with respect to the Fourier limit. 

We analyzed the values of F that we extracted from 
each fit used for Bloch frequency determination. A his- 
togram of the deviations from the Fourier-limited reso- 
nance Tp is reported in Fig. [9l Here one can observe 
a Gaussian-like distributed broadening of the AM reso- 
nance, with a mean value of about 1.5 mHz above the 
Fourier limit and a standard deviation of 1.3 mHz. If 
one considers a maximum line broadening of 2.8 mHz and 
assumes this estimate as the decoherence rate, it means 
that the coherence time is better than 330 s, which is 
again consistent with coUisional effect and with the BO's 
estimate as represented in FigUl 

Another source of sequential tunneling can be due to 
spurious AM . As shown in Sec. lIIIl the servo- loop system 
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controlling the lattice laser intensity has been optimized 
in order to avoid spurious amplitude modulation . Resid- 
ual amplitude noise on the lattice laser intensity induces 
resonant tunneling at the harmonics of the Bloch fre- 
quency. The level of the residual relative intensity noise 
is about -120 dB, which yields a fluctuation on the atomic 
size AfTij/AT ~ 10~^^m. This effects may limit the sen- 
sitivity on the Bloch frequency at the level of 2x 10~^ in 
relative units. 



C. Use as a force sensor: systematic effects 

We analyze in this section the systematic effects which 
affect the measurement of the gravity acceleration g. 
These effects are common to both BO measurement and 
WS resonant tunneling techniques. 

1. Lattice light shifts 



Agij and Agk- These sources of systematic shift must be 
estimated for a precise determination of the local gravity 
value g. We accomplished this task by a precise determi- 
nation of both the geometry of the incoming and reflected 
trapping beams and the position of the cloud with respect 
beam waist, with a relative uncertainty of 1%, as show in 
Fig. fTOl a). From this measurement we estimated a beam 
size on the atomic cloud w{zq) = 557(7) fj,m, considering 
also the beam ellipticity due to different beam sizes in 
the two transverse directions. 

An independent determination of the transverse beam 
size at the atomic sample position has been performed 
by measuring the radial atomic trap oscillation frequen- 
cies through parametric heating technique [5l|, 113 • On 
Fig. rTOT b) we recorded the radial characteristic frequency 
I'li — 5.3(1) Hz for an input laser power Pin = 0.848 W. 
Since the relation between these parameters is known 
[57| . we can express the potential depth Uq as function 
of the input power as 



The main contribution to systematic shift in local grav- 
ity measurement with trapped neutral atoms is due to 
the lattice light itself. The spatial inhomogeneity of both 
the intensity and the wave vector of the Gaussian lattice 
beam yields space-dependent variations of the potential 
energy, which results in 



Ua = kP„ 



and 



U{z) = UoI{z) (40) 



Utotiz) = Us{z) + Ui{z) cos(2fc(z)z) - mgz 



(37) 



where Us and Ui depend respectively on the squared sum 
and on the product of the electric field amplitudes of the 
two counter-propagating beams. Comparing this exter- 
nal potential with the ideal lattice in ([T|), besides the 
gravitational potential one recognizes another term due 
to the lattice intensity gradient. In the case of a red 
detuned lattice laser, atoms' wavefunctions are mainly 
located at the antinodes of the lattice where the inten- 
sity (and the potential) is maximum, so that the driving 
force is F = —dzV{z) = dz{Us{z) -\- Ui{z) — mgz). Fur- 
thermore, a change in the lattice constant d = 7r/fc(z) 
must be considered due to the Gouy effect [HI] , that is 



Hz) 



Mz)) 



(38) 



where $i(z) and ^r{z) are the phases for the incoming 
and reflected beams, respectively, and zq is the position 
of the atomic cloud. Thus, the gravity value becomes 



1 



2hkLU!B 



dUodz) 
dz 



hvB-^A<^Gouy{z) 



Z = Zq 

(39) 

where Uol{z) = Us{z) + Ui{z) is lattice field depth and 
A.<^Gouy{z) is the difference of the Gouy phase of the two 
beams. The Gaussian beam nature of the lattice intro- 
duces two extra terms in (|39p due to the intensity gra- 
dient and the Gouy phase shift that we call respectively 



where I{z) is the normalized intensity profile given by 
w and zq- For typical experimental parameters the two 
terms are Agu = 1.53(3) • lO'^m/s^ and Agu = 1.0(2) • 
10-**m/s2. 

It is possible to reduce systematic effects from the 
intensity gradient by employing a blue-detuned lattice 
laser [5a] . In this case, atoms are located at the inten- 
sity nodes where the light intensity is minimal so that 
the second term in (P7|) can be neglected. Since Us — U 
typically, the blue-detuned lattice light shift is reduced 
by a factor of two. 



2. Other sources of systematics and final uncertainty budget 

Beside the light shifts and the calibration of the lat- 
tice frequency, there are other technical and fundamental 
sources of shift of the Bloch frequency. 

As shown in the experimental apparatus description, 
during gravity measurement the vertical alignment of the 
lattice beam has been checked for each measurement. 
Moreover, during the measurement a tilt-meter with a 
resolution of 1.7 /irad and attached to the optical table 
has been employed to check alignment stability. This 
procedure has typical errors of 0.5 mrad leading to a rel- 
ative uncertainty 5g/ g — 10~® while the stability of the 
optical table over the duration of a typical measurement 
(< 1.5 hour) is better than 10 /irad. 

We estimated the contribution due to the distribution 
of masses around the experimental setup and we found 
that the biggest contribution to the local field is given 
by the 3.5 x 1.5 m optical table (20 cm thickness) that 
supports the experiment. Given the distance between 
the table surface of 25 cm and the weight of 640 kg we 
calculated the effect on the vertical direction to be = 
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FIG. 10. Characterization of the lattice beam propagation and intensity in the proximity of the atomic cloud, (a) The 
measurement of the beam width of two orthogonal transverse directions (respectively squares and circles) determines the 
trapping potential gradient (best fits are respectively the (blue) dot-dashed and the (red) dashed curves) , while (b) parametric 
excitation of radial motion is used to calibrate the intensity of the optical lattice at the place of the atoms. 



2 • 10"** m (2 . 10-9 in relative units). This effect 

is far below the current sensitivity. 

With the increase of the harmonic index £, it is pos- 
sible that for WS resonant tunneling measurements, at 
high modulation depths the contribution of off-resonant 
transitions has to be taken into account. The overlap of 
an infinite series of WS transitions induces a resonance 
asymmetry into the observed frequency spectrum and a 
consequent shift of the resonance peak. However the shift 
scales as {ujg t^)^^ and it is fully suppressed below 10~^° 
level. However, this effect couples with residual ampli- 
tude noise, which we already estimated to limit the sen- 
sitivity at the level of 20 ppb. Furthermore, off-resonant 
interband transitions can be non-negligible [58[ and with 
expected linewidths of the order of kHz. We experimen- 
tally checked possible shifts and broadening of the res- 
onance due to off resonant couplings and higher order 
effects by changing the modulation depth from 4% to 
10%. Within the statistical error in the determination of 
the resonance center and resonance linewidth we do not 
observe any shift or deviation from the Fourier limit in 
the measured linewidth. 

Another effect to be taken into account for these mea- 
surements comes from tidal forces. The peak-to-peak 
effect of tides at our site is of the order of 2 x 10~^m/s^, 
i.e., Aujb/uj < 100 ppb. This value is only two times 
smaller than the current best sensitivity for ujb measure- 
ment. Nevertheless the correction for this effect must be 
applied in accurate gravity measurements [l[ ■ Since each 
measurement lasts about 1 hour, the variation of g dur- 
ing every measurement due to tides is below 10~^m/s^ 
(i.e., below 10 ppb in relative units), then Earth tides 
does not affect the current uncertainty budget. 

In table |T] we present a summary of all the important 
systematic effects occurring in a gravity measurement are 



summarized. Any Bloch frequency determination is cur- 
rently limited by 200 ppb uncertainty due to systematic 
effects, mainly dominated by the lattice frequency uncer- 
tainty. Beside this effect, the force sensor based on WS 
resonant tunneling has proven a systematic error at level 
of 50 ppb. According to if a frequency stabilized 

lattice laser is employed, this level of sensitivity can be 
achieved after 26 s of modulation with the current setup 
at the 6-th harmonic. 



VI. CONCLUSIONS 

We explored the physical properties and the techni- 
cal problems behind the measurement of forces by atoms 
trapped in driven optical lattices. After fifteen years from 
the first detection [l3| of the quantum phenomenon of 
BO's, the precision and the accuracy of force measure- 
ment have reached the hundreds of ppb-level. This has 
been possible by exploiting the delocalization obtained 
by dynamical mechanisms and the matter-wave inter- 
ference nature of this phenomenon. In particular, we 
have shown that the DEBO measurement originates from 
modulation-induced delocalization leading to dispersion- 
less interference peak in momentum space [25| . 

Direct measurement of the Bloch frequency, and thus 
of the external forces, can be performed by means of in- 
situ measurements of the atomic cloud dynamics under 
AM or PM resonant tunneling among WS states. In par- 
ticular, our experimental and theoretical analysis shows 
that the AM resonant tunneling results the most accurate 
technique with higher potential for further improvement 
both in accuracy and in sensitivity. We provide a detailed 
characterization of the uncertainty budget of our force 
sensor. In the base setup, its uncertainty is due to the 
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Effect 


Correction (10 ^) 


Uncertainty (10 ^) 


Lattice wavelength fluctuations 





2 


Lattice beam vertical alignment 





0.1 


Inhom. Stark shift (beam geometry) 


14.0 -r-LI.O 


U.4 


Experiment timing 





0.2 


Tides 


-1.4 ^ 0.9 


<0.1 


Off-resonance tunneling 


<0.01 


0.2 


Systematics total 


17.2 ^ 22.5 


2.1 


(w/o lattice wavelength fluct.) 


II 


0.50 



TABLE I. Uncertainty budget concerning the gravity measurement with the atomic gravimeter. Correction values and their 
uncertainties depend on operating conditions. Typical values are given. 



lattice wavelength fluctuations. A frequency stabilized 
lattice laser would allow accurate measurements of forces 
with 50 ppb uncertainty. This value approaches per- 
formance demonstrated in free-fall atom gravimetry [H^ 
with the advantage of a small interferometer suited to 
microscale experiments. 
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Appendix A: Phase-modulation Hamiltonian 

In the rotating frame the Hamiltonian in Eq. ((T2|) be- 
comes 



TipM = Wrf'HpmZ^rf ^ ^^^RF^Z^RF ■ (Al) 



independent resonant terms 



The resulting expression 



-t-oo 



^PM =E l*">(*»l - E {^dcj 

71— — oo n— — oo 

X sin(a;j\/t — 0) exp [i/3£d/ (hujM) cos(a;jv/< — 0) — iwa/^] x 

,+f)(*„|+h.c.} (A2) 



+ 00 



->PM . 



can be further simplified by making use of the Jacobi- 
Anger expansion and by keeping then only the time- 



n— — on tr— — nn^ 



/2-0) 



\ Hum J \ ni^M J 



I +h.c. 



(A3) 



In order to retrieve the Hamiltonian in Eq. ([T^ one has 
to use the fact that Jo{x) + J2{x) ~ {2/x)Ji{x) and that 

^^--^\^\^-^"-{l/n.s for 'Jo 

with the coefficients eg defined in Eq. ([T5|) . 



Appendix B: Time-evolution operator 

It is convenient to express the WS states \'^n) in terms 
of the states defined in Eq. P7)) . which are eigenstates of 
the Hamiltonians in Eqs. ((TO)) and ((T3)) : 

= E e-""^^/'-'^) J-™(;^) \'^n-ram (Bl) 
m— — OO \ o / 

The time-evolved state at time t is given by 

oo 

= E (*™('5)l*n> e-""^'*/^ |<D,„(J)) = (B2) 

m— — OO 
oo 

^ ^-^nSt/l ^-^pin/2+4>) ^c,5t j^^^ | ^^^^^^ 

p,q— — oo 

with X — J'l/ [Fsd).. By using the addition theorem of 
Bcssel function, which states 

oo 

E e'''?"J,(x) J,+p(x) = e'P(--")/2 Jp[2 sin(«/2) x] , 

g— — OO 

(B3) 

one finds directly the expression in Ec^. ([1 
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Appendix C: Ballistic expansion at resonance 

Considering the special case on resonance, i5 = 0, the 
size a{t) of an atomic wave packet can be calculated in 
the asymptotic limit by means of the group velocities 
Vg(k), which arc defined, analogously to Eq. ([3]), as 



'.(fc) = ^^^^ = ^cos(£fcrf + 0), (CI) 



where -E'mod(fc) is the energy dispersion relation derived 
in Eq. (jl6p in the presence of resonant driving. Taking 
into account the uniform momentum distribution of the 
initial state |5'n), i-e., |(fc|4'„)p = 1, the variance at time 



t results 

aljt) fa J dx x"^ J dk S[x ~ Vg{k)t\^ = v\t^ 

(C2) 

in agreement with Eq. (^0)) . Since the computation of 
a\ {t) is rather lengthy, wc shall report only the final 
result for the AM modulation, 

al^ (t) = <2 smc[S/TF{t)f + 2V2 vit d ei/{hujB)-x 

X cos [0 - (wM + iuJB)t/2] sine [6/TF{t)] + (t% (0), 

" (C3) 

with the coefficients ei defined in Eq. ([15]), Tpit) the 
resonance linewidth, and CTv]/^ (0) the size of the initial 
state |^'„). 
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